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Abstract 

DAEDALUS, the Decay At-rest Experiment for Sep at a 
Laboratory for Underground Science will look for evidence 
of CP- violation in the neutrino sector, an ingredient in the- 
ories that seek to explain the matter/antimatter asymmetry 
in our universe. It will make a precision measurement of 
the oscillations of muon antineutrinos to electron antineu- 
trinos using multiple neutrino sources created by low-cost 
compact cyclotrons. The experiment utilizes decay-at-rest 
neutrino beams produced by 800 MeV protons impinging 
a beam target of graphite and copper. Two well established 
Monte Carlo codes, MARS and GEANT4, have been used 
to optimise the design and the performance of the target. A 
study of the results obtained with these two codes is pre- 
sented in this paper. 

INTRODUCTION 

Neutrino oscillation, the quantum mechanical process 
through which a neutrino created in one of three flavor 
states can be measured in another, is governed by three 
mixing angles, a CP-violating phase (Sep), the neutrino 
mass differences, and the ratio of the experimental baseline 
to the neutrino energy (L/E). All of these parameters, ex- 
cept Sep, have been measured. If Sep is non-zero, there 
is an asymmetry between how neutrinos and antineutrinos 
behave and this may help explain why there is more matter 
than antimatter in the universe. 

DAE<5ALUS [ 1 ] will search for evidence of CP-violation 
in the neutrino sector by measuring the oscillation of muon 
antineutrinos into electron antineutrinos. DAE<5ALUS will 
measure the oscillation probability at three distances and 
use these measurements to fit the probability formula and 
extract the 8c p parameter. The experiment will use three 
neutrino sources and one large underground detector to 
make these measurements. Each neutrino source will con- 
sist of a compact cyclotron that will output 800 MeV 
protons into a beam target made of graphite and copper. 
While the graphite is the primary producer of pions in our 
target design, the copper effectively reduces backgrounds 
and transfers heat to keep the graphite below melting tem- 
peratures. 

To design the DAEJALUS targets, we simulate and com- 
pare the neutrino flux produced by various beam target ge- 
ometries. Because the neutrino production rate, thus the 
pion production rate is our figure of merit, we must en- 
sure that our simulations agree with experimental measure- 
ments of pion production. In this paper, we present a com- 



parison of 7r + production in the simulation codes MARS 
and GEANT4 to available experimental data for protons 
with energies between 300 and 2500 MeV impinging on a 
graphite target. 

TARGET DESIGN 

Proton collisions in the carbon target will create A res- 
onances which then decay to create pions. The tt + decay 
to create muon antineutrinos, the source of the oscillation 
signal, and the 7r~ can decay to create electron antineu- 
trinos, the dominant background for the experiment. Pi- 
ons produced in the graphite travel into the copper prior 
to decaying. Once in the copper, the ir~s can be cap- 
tured, preventing their decay. This reduces the ratio of elec- 
tron antineutrino background to muon antineutrino signal 
to - 10- 3 0. 

There are many interactions that occur in the beam target 
which ultimately affect the production of neutrinos. Pro- 
tons lose energy as they traverse the target material prior to 
interacting, potentially preventing them from producing A 
resonances. Protons of sufficient energy can also re-interact 
causing multiple A resonances. Before considering these 
effects, it is important to benchmark the cross-sections that 
MARS and GEANT4 use for pion production in proton on 
carbon (p — > C) collisions. 

SIMULATION 

In experimental studies, p —> C inelastic cross-sections 
and pion production cross-sections are measured using a 
thin carbon target, with thickness between 1-1.79 g/cm 2 
or 4.5 - 8 mm Q. We simulate a 5 mm target 
for comparison to the experimental literature. While the 
DAE<5ALUS targets will be on the order of 1 m thick, when 
benchmarking the two codes with a target thickness much 
smaller than the interaction length of the incident parti- 
cles, we can neglect effects such as energy loss and re- 
interaction in the target. 

In both the MARS and GEANT4 simulations, we count 
the total number of ix + produced and number of inelastic 
proton interactions. We use these values and the simulation 
geometry and proton beam parameters to derive the cross- 
sections we compare to experimental data. 

MARS 

MARS |6| uses the Cascade Exciton Model (CEM) for 
materials with atomic mass A > 3 and energy less than 
5.0 GeV. CEM considers nuclear reactions as proceeding 



through three states - cascade, pre-equilibrium, and equi- 
librium |7|. First, there is the intra-nuclear cascade where 
cascade particles are emitted leaving the residual nucleus in 
an excited state. Second, the nucleus relaxes according to 
the exciton model of the pre-equilibrium decay. Third, the 
equilibrium evaporative stage of the reaction occurs [7 1 [8 1. 

GEANT4 

Three independent models have been used in 
GEANT4 [9 1 to simulate the inelastic proton interac- 
tions: the Bertini model, the Binary Cascade and the Liege 
(INCL) intra-nuclear cascade model coupled with the 
independent evaporation/fission code ABLA. 

The Bertini Cascade Model [ 10 1 generates the final state 
for hadron inelastic scattering by simulating the intra- 
nuclear cascade. Incident hadrons collide with protons 
and neutrons in the target nucleus and produce secondaries 
which in turn collide with other nucleons, the whole cas- 
cade being stopped when all the particles which can escape 
the nucleus have done so. Relativistic kinematics is ap- 
plied throughout the cascade and the Pauli exclusion prin- 
ciple and conformity with the energy conservation law is 
checked. This model has been validated up to 10 GeV in- 
cident energy and is performing well for incident protons, 
neutrons, pions, photons and nuclear isotopes. 

In the Binary Cascade Model ifTTl the propagation 
through the nucleus of the incident hadron and the secon- 
daries it produces is modeled by a cascade series of two- 
particle collision, hence the name binary cascade. Between 
collisions the hadrons are transported in the field of the nu- 
cleus by Runge-Kutta method. The model is valid for inci- 
dent protons, neutrons and pions and it reproduces detailed 
proton and neutron cross section data in the region 0-10 
GeV. 

To respond to the increasing user requirements from 
the nuclear physics community, the GEANT4 collabora- 
tion set a goal to complement the theory-driven models in 
this regime (the Bertini cascade and Binary cascade being 
the most widely used) with the inclusion of the INCL code 
also known as Liege cascade, often used with the evapo- 
ration/fission code ABLA fl2l . The code was validated 
recently against spallation data. It supports projectiles like 
protons in the energy range 200 MeV - 3 GeV. 

RESULTS 

We provide two comparisons of simulation results to ex- 
perimental data. First, we compare the p — >• C total inelas- 
tic cross-sections, crinelastio and the tt + production total 
cross-sections, a n + , for GEANT4, MARS and experiment. 
Second, we compare the ratio of these cross-sections. In 
this paper, we compare to the experimental pion production 
cross-sections from references fl3], |@], iTOl . lfj"4"l . and [1T51 
and the experimental p — !> C inelastic scattering cross- 
section from references [16] and ifPTl . 

For the proton inelastic cross-section, the three theo- 
retical models used in GEANT4 predict the same values. 



These values are compared with the MARS model and ex- 
perimental data in Fig. Q] The GEANT4 models closely 
align with data for energies less than 1 GeV. Above this 
energy, the data is highly varied, but is centered around 
both the GEANT4 and MARS models. Since these varying 
features do not reflect nature but rather experimental error, 
we use the GEANT4 models' prediction to represent the 
experimental data in the second half of our results. 
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Figure 1: Plot of experimental and theoretical models' pre- 
dictions for CTinclastic- 

Out of the three GEANT4 models investigated, the 
Bertini model predictions for the total pion production 
cross-section are closest to the available experimental data. 
The comparison of the MARS and GEANT4 model predic- 
tions with experimental data is shown in Fig. [2] Both codes 
compare well to the experimental values. We fit a third or- 
der polynomial to the experimental data to use in our final 
cross-section comparison. 
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Figure 2: Plot of experimental and theoretical models pre- 
dictions for ov+ . 

The ratio of the number of protons producing 7r + and 
the number of protons which interact inelastically in our 
thin target is given by: 
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where n is the number density of carbon atoms and x is 
the thickness of the target. To obtain the righthand ap- 
proximation, we Taylor expand around x over l/^ov) and 
l/(na p ), the interaction lengths of the two respective pro- 
cesses, to second order. We can make this approximation 
because the target thickness (5 mm) is much smaller than 
the interaction length of either process (^1 m). 

In Fig. [3] the GEANT4 and MARS cross-section ratios 
are obtained by dividing the simulation values in Fig.|2]by 
the simulation values in Fig.Q] The "Ratio of Fits to Data" 
is obtained by dividing the third order polynomial fit to the 
experimental data in Fig. |2]by the GEANT4 model predic- 
tion in Fig.[T] 
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Figure 3: Plot of MARS and GEANT simulation results 
and the ratio given by the two fits through the experimental 
data. 



CONCLUSION 

Pion production calibration is crucial for ensuring ac- 
curate simulation of the neutrino fluxes created by the 
DAE<5ALUS sources. Calibrating for energies lower than 
the initial proton energy is important for simulating the 
pion production in thicker targets where more protons in- 
teract but do so at lower than their initial energies. 

From this study we find that MARS and GEANT4, 
for the models presented, consistently agree on pion pro- 
duction rates within 30% and that both codes compare 
well to experimental data in the energy region of interest, 
< 800 MeV. The information from this study can be used 
to reweight simulation data, such that simulated flux more 
accurately reflects experimental data. 

The neutrino flux predictions we make with these sim- 
ulations will be used to decide on experimental param- 
eters such as cyclotron power cycle and determine the 



required running-time needed to achieve DAEi5ALUS's 
physics goals. 
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